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SUMMARY 

The nonlinear discrete-vortex method has been extended to treat the problem of 
asymmetric flows past a wing with leading-edge separation, including steady and unsteady 
flows. The problem is formulated in terms of a body-fixed frame of reference and the 
nonlinear discrete-vortex method is modified accordingly. Although the method is general, 
only examples of flows past delta wings are presented due to the availability of experi- 
mental data as well as approximate theories. Comparison of our results wi+h the experi- 
mental results of Harvey for a delta wing undergoing a steady rolling motion at zero 
angle of attack demonstrates the superiority of the present method over existing appro- 
ximate theories in obtaining highly accurate loads. Numerical results for yawed wings 
at large angles of attack are also presented, In all cases, total-load coefficients, 
pressure distributions and shapes of the free-vortex sheets are shown. 
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LIST OF SYMBOLS 


AR aspect ratio 

b(x) local half span 

C rol ling -moment coefficient 

A 

Cp pressure coefficient 

AC net pressure coefficient (C - C ) 

fJ pi P2 

Cj, root chord 

ACj, winglet root chord, a characteristic length 

e" unit vector in the direction of the freestream velocity 

oo 

T,T,¥ unit vectors of wing-fixed frame of reference 

i length of a vortex segment 

£ total number of vortex segments of the model 

nj^ unit normal to the wing surface 

"n, unit normal to the v/ake surface 

w 

7 . position vector of a field point 

S^7) wing surface 

;t dimensionless time 

U freestream velocity 

w{"r,t) wake surface 

x,y,z wing-fixed frame of reference 

a, a angle of attack and rate of pitch 

angle of yaw and rate of yaw 

Y,Y angle of roll and rate of roll 

w angular velocity of wing 

n frequency of rolling velocity 

<j) disturbance velocity potential 


Subscrii 


i source point 

j field point 

k time-step number 

LE leading edge 
TE trailing edge 
t tangent 
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I. INTRODUCTION 


In recent years, development of analytical and numerical methods for predicting 
the aerodynamic characteristics of wings exhibiting leading-edge and/or wing-tip 
separation has received considerable attention. The literature contains several 
analytical methods which are based on simplifying assumptions. These include the assumption 
that the flow is conical and the assumption that the axial gradients are much smaller 
than the lateral gradients of the flow properties {slender-body theory). These assump- 
tions represent inaccurate modelling of the full three-dimensional flow and violate the 
trail inq-edge Kutta condition. Methods of this type were presented by Brown and Michael S 
Mangier and Smiths and Smith^. 

With the advent of high speed computers, new techniques which avoid these simplify- 
ing assumptions were developed. The nonlinear, discrete-vortex technique is among 
those successful techniques capable of calculating full three-dimensional flows. With 
this technique, the exact governing equation and the corresponding boundary conditions 
are satisfied and hence the solution of the inviscid problem is exact in that sense. It 
is not restricted by the shape of the wing planform or the range of angle of attack as 
long as vortex breakdown does not occur in the vicinity of the wing (stall phenomenon of 
wings of low-aspect ratio). However, the separation line is assumed to be known a 
priori along the sharp edges of the wing. Methods of this type were developed for 
steady symmetric flows by Belotserkovskii^, Rehbach®’®, Kandil’, Kandil, Mook and 
Nayfeh®"“ . 

Due to the accuracy and simplicity of the nonlinear discrete- vortex technique, it 
was extended to treat unsteady symmetric flows by Belotserkovskii and Nisht^^, Atta^^, 

Atta, Kandil, Mook and Nayfeh^ ^ ® . The source of unsteadiness in the flow may be 
general, e.g. a sudden translational acceleration of the wing^^, a gust wind which 
changes the wing angle of attack* > or an oscillatory pitching motion of the wing*®. 

The method is characterized by its capability of obtaining the transient as well as the 
steady-state aerodynamic characteristics of the wing. The technique was also extended 
by Kandil, Mook, and Nayfeh*® to treat steady asymmetric flows past a large aspect-ratio 
rectangular wing. This case is a simulation of the problem of aerodynamic interference 
which arises when a small aircraft penetrates the wake of a large aircraft. Although 
the authors didn't account for the wing-tip separation of the trailing wing (due to its 
large aspect ratio), they did account for the wing-tip separation of the leading wing. 

In the present paper, this technique is applied to steady and unsteady asymmetric 
flows past highly swept-back wings with sharp-edges. Delta wings are chosen as numeri- 
cal examples due to the availability of experimental and theoretical results. Another 
reason for this choice is that it represents a severe numerical test of the technique 
ov/ing to the presence of the vortical spiral cones which emanate from the leading edges 
and extend over a large portion of the wing surface. Moreover, the free-vortex surfaces 
are represented by a series of segmented vortex lines which approach the bound- vortex 
lattice representing the wing surface during the development of the numerical solution. 
Hence, strong singularities may arise due to the interaction of close vortex lines and 
therefore safeguards must be imposed to eliminate such singularities. 

The problem of steady asymmetric flows was treated earlier by Pullin*^, Hanin and 
Mishne*®, Jones*®, and Cohen and Nimri^“. Specifically, the steady flow past a yawed 
slender delta wing ms considered in references 17 and 19 while the flow past a slender 
delta wing rolling steadily was considered in references 18 and 20. These theories are 
based primarily on the approximations of slender-body theory. The method of Brown and 
Michael was extended for the rolling wing*®’^® while the method of Mangier and Smith and 
the improved method of Smith were extended for the yawed wing*^’*®. 
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Therefore, with these approxim?.tions as well as the simplified approximate modell- 
ing of the vortex sheets shed from the leading edges, one may expect substantial differ- 
ences between the predicted and the experimental results. In fact, in all asymmetric 
results of these theories substantial errors exist in predicting the suction peak of the 
pressure near the leading edges. This was the case on the receding face of a rolling 
delta wing at zero angle of attack as reported in references 18 and 20,. For the case of 
a yawed wing at high angle of attack, large errors were reported in predicting the 
suction peak of the pressure on the windward side of the wing^'^’^®. In the latter case 
the errors were attributed to the substantial secondary separation which develops due to 
the adverse pressure gradients on the windward side of the wing. 

From our point of view, this is only one portion of the cause because as it is 
well-known the problem of secondary separation is a viscous phenomenon and it cannot be 
treated by an inviscid model. The other portion of the cause is in fact the modelling 
of the separated flow and the slender-body assumption . Moreover, on the leeward side 
of the wing, the core of the primary vortex moves outboard as the angle of yaw is in- 
creased. Hence the suction peak on this side disappears and so does the adverse pres- 
sure gradients. Thus, secondary separation diminishes and one can expect an inviscid 
three-dimensional model to yield highly accurate results on the leeward side. 

On the experimental side, Fink^^ and Harvey^^ considered steady flows over yawed 
slender-delta wings. Later, Harvey^^ considered flows over a steadily rolling delta 
wing. Pressure distributions, local rolling-moment coefficients, local normal -force 
coefficients, total -load coefficients and positions of the vortex cores were reported in 
these experiments. These data are used by many investigators for checking the accuracy 
of their theories. In the present paper, we also consider the same data to check our 
results for steady asymmetric flows. We also compare our results with available approx- 
imate theories. The method is also extended to the problem of unsteady asymmetric flows 
past a rolling wing. 
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II. FORMULATION OF THE PROBLEM 


We consider a thin delta wing in a uniform stream and let U be the free stream 

DO 

velocity and oxyz be a wing-fixed frame of reference. The wing edges lie in the xy- 
plane, the x-axis bisects its apex angle and the xz-plane is its plane of symmetry. 

Euler's angles a, 3 and y are used to define the angle of attack, the angle of yaw and 
the angle of roll of the wing, respectivelyj^ see Figure 1. To construct these angles in 
this order, we start from a position where is parallel to the x-axis and successively 

allow for the_ positive rotations a,S,y about the y,z, and x axes, respectively. The 
unit vector in the direction of the freestream velocity is expressed in terms of 

these angles and the base unit vectors o*' the fixed- frame of reference by 

e^ = cosa cos3i + (sina siny - cosa sinB cosy)j + fsina cosy + cosa sinB siny)l< (1) 

Next, we assume that the wing is rotating with an angular velocity cj which can be express- 
ed in terms of Euler's angles and their rates of change cs 

Hi = w^i + Wyj + to^k 

= (a sin3 + y)T + (a cos3 cosy + 3 slnr)f + (-a »- 0 Sp si"y + 3 cosy)k (2) 

The fluid flow is assumed to be ideal. The assimption that the wing edges are 
sharp fixes the separation lines along these edges. Vorticity is shed from these edges 
in the form of free surfaces of tangential discoptinuiiy {free-vortex sheets). Moreover, 
the flow is assumed to be irrotational in the region R exterior to the wing and its 
free-vortex sheets. Accordingly, the flow in R is governed by Laplace's equation 

V^(i) = 0 (3) 

where 4){F,t} is the disturbance potential. On the boundary 3R, tf) satisfies the following 
boundary conditions. The flow must be tangent to the wing surface; that is, 

+ V(j) - u X r) • = 0 on S(r) = 0 (4) 

where = VS/|VS| is the unit normal to the wing surface S. The pressure is continuous 

and no flow exists across the free-vortex sheets; these dynamic and kinematic conditions 
yield 

" Sr Sa " ~ ^ - 2 = 0 

on w(F,t) = 0 

T ^ I f 

(F^+ V<|) - w X r ) " F^ - a t ° w(r,t) = 0 

where the subscripts 1 and_2 refer to the upper and lower surfaces of the free-vortex 

sheets, respectively, and n,, = Vw/IVw] is the unit normal to the free-vortex sheets. 

w 

The Kutta condition must be satisfied along the edges of separation, that is the flow 
leaves smoothly off the edges and the pressure is continuous across these edges. These 
are expressed by 

(F^ + V(|) “ to 'P') • % = 0 on S! =0 

TEjLE 


(5) 

( 6 ) 


(7) 
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and 


ACp = {V(f)i - V<i) 2 ) • [2(w X r - e^) - V(f)i - V(}) 2 ] - 2 (<j>i - <j) 2 ) = 0 

on S[ =0 (8) 

TE.LE 

Far from the wing and its free-vortex sheets, the disturbance velocity vanishes; that is, 

V(j) -> 0 (9) 

The complexity of the problem stems from the boundary £onditions (5) and (6). These 
conditions are to be satisfied at the_free“Vortex sheets w(r,t) which are unknowns of the 
problem. In fact, both w(r,t) and <{)(r,t) are dependent upon each other and hence a 
numerical solution is fruitful for this situation. In the next section, we discuss the 
method of sulution of the steady problem and describe its extension to the unsteady 
problem. 
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III. METHOD OF SOLUTION 


1 . Nonlinear Discrete-Vortex Method for the Steady Problem 

The method of solution of the steady asymmetric flow is obtained by generalizing 
the approach of Kandil’ and Kandil, Mook and Nayfeh® for the steady symmetric flow. For 
the sake of completenessj we outline the basic approach and then describe the generalization 
to the current problem. 

In the basic approach of the lifting problem, a lifting surface may be replaced by 
a bound-vortex sheet with unknown strength. The velocity potential of this sheet satis- 
fies equation (3) in R and the boundary condition (9). Furthermore, this continuous 
vortex sheet can be accurately approximated by a lattice of bound-vortex filaments with 
unknown circulations (bound-vortex lattice) provided it is constructed within certain 
rules. This point will be undertaken later in detail in the next section. 

Along the edges of separation no bound-vortex segments are placed; otherwise the 
Kutta condition, equation (8), will be partially violated. At these edges, the starting 
vortex is shed and convected with the local velocity and in a steady flow no more vor- 
tices are shed thereafter. According to the theorem of spatial conservation of circu- 
lation (Kelvin's theorem), the ends of the bound-vortex lattice lines closest to the 
edges are connected to vortex lines which extend downstream to infinity where the start- 
ing vortex is assumed to be. JThese vortex lines are called free-vortex lines and repre- 
sent the free-vortex sheet w(r). So far, the model satisfies equations (3) and (9) and 
partially satisfies the Kutta condition. 

For steady symmetric flows, equations (4) - (8) immediately yield the corresponding 
boundary conditions upon setting w = 0 and dropping the time-dependent terms while' 
equation (1) yields the corresponding free stream velocity upon setting (3 = ■y = 0. With 
the discrete metjiod, the resulting equations are sat^isfied at certain points on the 
known surface S(r) and the still unknown surface w(r). This is achieved by successive 
iterative cycles. In the first st£p of the cycle, we satisfy the flow tangency condition 
on S(r) with an assumed surface w(r) to find a circulation distribution T. In the next 
step, we satisfy the kinematic and dynamic boundary conditions on w(r) by using the r_ 
distribution and find w(r). These cycles are repeated until the r distribution or w{r) 
do not change within certain prescribed tolerances. The Kutta condition is then satis- 
fied automatically. 

For steady asymmetric flows, equations (4) - (8) also yield the corresponding 
boundary conditions upon dropping the time dependent terms. Now, to obtain a steady flow 
it is_a necessary condition, although not sufficient, that the angular velocity of the 
wing m be uniform. Moreover_^ the magnitude and direction of w as well as the orientation 
of the wing vnth respect to e^ (as defined by a,3,Y) niust have values such that no time- 

dependent disturbance is generated in the flow. For a thin flat wing with m = 0 (a = 

S = Y = 0), nontrivial, steady, asymmetric flows occur in three cases. First, y = 0 and 
the free stream velocity is given by 

e"^ = cosa cosS T - cosa sinB T + sina 1< (10a) 

This case represents the steady flow past a yawed wing at an angle of attack. It was 
considered in references 17, 19, and 22. Second, 8=0 and the free stream velocity is 
given by 

e^^ = cosa T + sina siny j + sina cosy k (10b) 
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This case represents a steady flow past a wing at a banking angle and an angle of attack. 
Third, a = 0 and the free stream velocity is given by 

■e^ = COS3 T - sing cosy J + sing siny 1< (10c) 

This case represents a steady flow past a yawed wing at a banking angle. Fourth, when 
all the angles a, 3, and y are different from zero and the free stream velocity is given 
by egi.ation (1). All the cases considered above can be treated by the vortex lattice 
method as reported in references 7-10. 

whi the angular velocity w is not equal t£ zero, there are still two cases where 
the flow is 3 :eady. First, a = 3 = « = 3 = 0, w = yi and the free stream velocity is 
giv^'n oy 

'i, "1 (lla) 

This case represents a steadily rolling wing about its x-axis at a. zero angle, of at;tack. 
It, was considered in references 18, 20, and 23. Second, B = y = a = 0, w = yi + 3k, tana 
= 3/y and the free stream velocity is given by 

= cosa T + sina k" (11b) 

This case represents a steadily rolling wing about the wind axis at an angle of attack. 

We can easily see that the first case is a special case of the present one. This case was 
considered in reference 20. 

When applying the steady version of the discrete-vortex method to the last two 
cases, one must use a wing-fixed frame of reference; otherwise, the flow would no longer 
be steady. 

2. Nonlinear Discrete-Vortex Method for the Unsteady-Flow Problem 

The method of solution for unsteady asymmetric flows is obtained by generalizing the 
approach of Atta^^ and Atta, Kandil, Mook and Nayfeh^'^*^^ for unsteady symmetric flows. 

In these references, a space-fixed frame of reference V'/as used, the wing was taken to be 
fixed in the flow, and the source of unsteadiness was introduced through the free stream 
velocity. In the present paper, a wing-fixed frame of reference is used, the wing is 
rotating at a nonuniform angular velocity w(t), and the free stream velocity is uniform. 

In either case, the bound circulation around the wing continuously changes and this 
is accompanied by a continuous process of formation and shedding of vortexes from the 
edges of separation to restore the smoothness of the flow at the edges (Kutta condition). 
Within any infinitesimal time step, the change in the bound circulation around the wing 
is met by the formation of an infinitesimal vortex strip emanating from an edge of 
separation which has a strength of equal and opposite sense to the change of the bound 
circulation. This shed vortex is convected downstream with the local particle velocity. 
Hence a vortex sheet is continuously growing dov/nstream as long as the unsteadiness of 
the flow prevails. 

Now if the continuous motion of the wing is discretized into a series of impulsive 
changes occuring at discrete time steps, the continuously growing vortex sheet can be 
replaced by a growing vortex lattice in the viake. This is the salient difference between 
the unsteady and steady flow models. At each time step, we solve the problem using a 
method similar to that of the steady flow. Here, the boundary conditions, equations 
(4) - (8), must be satisfied at each time step. In the next section, we consider in 
detail the implementation of this method. 
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IV. IMPLEMENTATION OF THE METHOD 

1 . Construction of the Discrete-Vortex Model 

In Figure 2, we show how the discrete-vortex model is constructed for a delta wing. 
Although the example discussed here is for a thin, flat, delta wing, the method is 
general and is not restricted by the geometrical parameters of the wing; e.g. camber, 
aspect or thickness ratios or wing planform. 

The first step is to divide the wing into rectangular and cropped-delta winglets as 
shown by the dashed lines in Figure 2. a. A rectangular winglet is aerodynamically 
represented by a spanwise bound- vortex segment of constant circulation T^. . This segment 

is placed at the quarter-chord length of the winglet (the chord length of the rectangular 
winglet is the characteristic length of the problem). In addition a control point is 
placed at the three quarter-chord length. The choice of these positions is suggested by 
thin airfoil theory^'*. It can be shown that the bound-vortex sheet representing the two- 
dimensional flow around a flat plate at an angle of attack can be replaced by a point 
vortex of the same strength as that of the continuous vortex sheet under the following 
conditions: a) the point vortex is placed at the quarter-chord length and b) the flow 

tangency condition is enforced at only one point at the three-quarter-chord length. 

On the other hand, a cropped-delta winglet is aerodynamical ly represented by a 
bound-vortex segment of constant circulation. This vortex segment is directed along the 
perpendicular from the midpoint of the winglet root chord to its leading edge. With this 
choice it can be seen that the vorticity of this vortex segment does not have a component 
along the leading edge and hence the Kutta condition is approximately satisfied along 
this edge, 

Chordwise bound-vortex segments arise due to the differences in the strengths of the 
neighboring spanwise, bound-vortex segments. In this way, a bound-vortex lattice which 
replaces the continuous, bound-vortex sheet is constructed. The model is completed by 
adding free-vortex lines, representing the continuous free-vortex sheets at the ends of 
the bound-vortex lattice along the edges of separation - the leading and trailing edges. 
Each line is divided into a series of small, straight segments (near-wake region) and one 
semi-infinite vortex line (far-wake region). The upstream end of each segment represents 
a control point of the wake surface where the kinematic and dynamic boundary conditions 
are satisfied. The initial positions and shapes of these lines are prescribed. The 
resulting model is shown in Figure 2b. This model has an unknown circulation distribu- 
tion and a wake that can be deformed to satisfy the boundary conditions. 

The model described above is used to solve the steady-flow problem by satisfying the 
corresponding boundary conditions given in Section III.l. On the other hand, if the 
problem under consideration is for an unsteady flow which starts from a steady flow 
situation, then the solution of the steady-flow problem serves as an initial condition to 
the unsteady problem. Furthermore, if the problem under consideration is for an unsteady 
flow which arises from an impulsive motion of the wing, then the initial condition 
corresponds also to the solution of the model given above, but with the wakes removed 
from the model . 

2. Calculation of the Velocity Field 

To satisfy the boundary conditions on the wing and its wake and to calculate the 
surface pressure distribution, one needs an accurate method to calculate the velocity at 
any field point "r. at any time step If the field point is off the wing and its wake, 

then the velocity is given by 
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(12) 


V(rj.tk) = X r. 

where 


is the induced velocity from all the vortex segments of the model. The parameters on the 
right-hand side of equation (13) are those of Biot-Savart' s law^'*. The number of vortex 
segments n(t|^J is a function of the time step tj^ due to the growing vortex lattice in the 
wake in the unsteady-flow problem. To avoid extremely large induced velocities, an 
"artificial viscosity" is introduced (as in ref. 25) in the form of an exponential 
multiplier which causes the induced velocity to approach zero as the vortex is approached. 

When the field point is on the wing surface or on the wake surface, one has to 
account for the induced tangential velocity due to the local strength of the vortex 

sheet. In Figure 2.c, we show the parameters involved in calculating the components of 

the induced tangential velocity in the x and y directions at a point p for a quadrila- 
teral vortex element in the xy-plane. With linear interpolation, it is easy to show that 
these components are given by 

\^{z = 0 ") = ± 0/2 Jl?)[ri{xi + - x) + r,(x - X:)]ei (14a) 

\y{z = Q-) = + (1/2 £|)[rn(yi + U - y) + rs(y - yi)}e2 (14b) 

- 

where the arguments z = 0 and z = 0 correspond to the upper surface and the lower 
surface of the wing, respectively. Equations (14) must be added to equation (12) if one 
is to calculate the pressure distribution on the upper and lower surfaces by using 
Bernoulli's equation. Extension of equations (14) to a general, quadrilateral vortex 
element is straightforward. 

3. Implementation of the Boundary Conditions 
a. Steady-Flow Problem 

The boundary conditions on the wing surface S(7) and the wake surface w(F) are 
satisfied by an iterative process. To initiate the iterative process, one needs to 
prescribe an initial geometry of the wake surface. It has been found from several 
numerical tests that the number of iterative cycles required to achieve the solution can 
be reduced by an appropriate choice of the initial geometry. This initial geometry 
depends on the problem under consideration and thus it varies from one problem to the 
other. 

For instance, the number of iterative cycles for the steady, symmetric-flow problem 
is reduced by about 20% when the free-vortex lines emanating from the leading edge are 
prescribed to be straight lines pitched at one half the wing angle of attack. In addi- 
tion, those lines emanating from the trailing edge are straight lines pitched at one 
third the wing angle of attack. Here, the comparison is made with respect to the number 
of iterative cycles required for the same problem when all the free-vortex lines are 
prescribed to be also straight lines but are pitched at an angle equal to the wing angle 
of attack. 

In the case of a steadily, rolling wing at zero angle of attack, an appropriate 
initial guess is found to be related to an angle 0(r) = + ^ tan~M«xrl/U . Here, we 




= I [ri(\)/4TTh.j(rj.,t,^)][cos0,.(r^,t,^) . cosOz . (r^. ,t^)]e^ (r - ,t^) (13) 
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specify the free-vortex lines emanating from the edges of the advancing side to be 
straight lines pitched at the angle + B while those emanating from the edges of the 
receding side to be straight lines pitched at the angle - 6. 


Next, the flow-tangency condition and the spatial conservation of circulation are 
applied at the control points and node points, respectively, of the bound-vortex lattice. 
Thus, m obtain a set of linear algebric equations” which yields the circulation dis- 
tribution V.. 

With the circulation distribution known, the kinematic and dynamic boundary condi- 
tions at the control points of the free-vortex lines are satisfied. For steady flows, 
these two conditions are combined into a simple condition in which we require that each 
vortex segment in the wake be aligned with the local velocity at its upstream end (a 
control point on the wake surface). This means that each vortex segment is a segment of 
a streamline (kinematic condition). Moreover, it means that the force on each vortex 
segment is zero according to Kutta-Joukowski theorem in the small (dynamic condition). 
This process is carried out by calculating the downstream end of each vortex segment 
according to 


'"j+1 


= '■j " 


(15) 


where and are the position vectors of the upstream and downstream ends, re- 
spectively, is the segment length and is the velocity at its upstream point 
[equation (12) for steady flows]. 


The iteration scheme moves back and forth from the control points of the bound- 
vortex lattice to the control points of the free-vortex lines until convergence is 
achieved. We consider the iteration scheme converged when the variation in the circu- 
lation distribution or the displacement of the downstream ends of the free-vortex segment 
between two successive iteration cycles does not exceed a certain prescribed tolerance. 
Once convergence is achieved, we calculate the pressure distribution and the total load 
coefficients. 


b. Unsteady-Flow Problem 

Here, we consider the problem of unsteady flow which starts from a steady flow 
situation. We recall from Section III. 2 that the continuous motion of the wing is 
discretized into a series of impulsive changes occuring at discrete time steps. At each 
time step t|^, a set of starting vortices develops along the edges of separation and are 

shed with the local velocities to restore the smoothness of flow at the edges (Kutta 
condition). In the same time, the starting vortices shed in the wake at earlier time 
steps are convected downstream with the local velocities without changing their strengths. 
This process satisfies the kinematic condition on the wake (a wake element moves along 
the direction of the local velocity) and it also satisfies the dynamic condition on the 
wake (a wake element satisfies Kelvin-Helmholtz theorem). 

The position of any shed vortex 7^. at any time step tj^ is determined by 

Vj(tk) = + (tk - V,) (16) 

where t|^_-| is the preceding time step and is given by equation (12). The strength of 

any newly shed vortex is related to the change in the bound circulation. Hence, with the 
positions of the shed vortices known from equation (16) and with the strength of the 
newly shed vortices given in terms of the change of the bound circulation, the flow 
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tangency condition at the control points of the wing yields the unknown circulation 
distribution. To account for the error in equation (16) because of using the velocity 
V,j(7j ,tj^_^ ) at the preceding time step ^ rather than the current time step t^^, an 

iteration procedure similar to that of the steady-flow problem is performed. In this 
regard, an alternative equation was given by Summa^'"’. 

In both the steady and unsteady flows, the only difference between the symmetric and 
asymmetric problems is the longer computational time required for the latter problem as 
compared to that of the former problem. In the former problem, we need only to use half 
the wing to obtain the solution because of the symmetry of the flow. In the latter 
problem, the whole wing must be used to obtain the solution. 

4. Calculation of the Pressure Coefficients 


The distribution of the pressure coefficient on the upper and lower surfaces of the 
wing is calculated by using Bernoulli's equation in terms of a wing-fixed frame of 
reference^ 


Cp{r^^t[^) = - [v(rt,t|^)]2 + 2vTr^.,t|^) • [m(t|^) x r^ - ej - 2 




3t 


(17) 


where 


r"'^ is the position vector of the control point, the positive and negative superscripts 
« — * 
refer to the upper surface and lower surface of the wing, respectively, and VtJ), V 

T A 

and are given by equations (13), (14a) and (14b), respectively. The pressure is 

calculated at the control points of the bound-vortex lattice because these are the 
points where the flow tangency condition is enforced. 


In the steady-flow problem, the last term of the right-hand side of equation (17) 
is zero and all the other terms are time independent. Hence, calculating the pressure on 
the upper and lower surfaces is straightforward. But in the unsteady-flow problem, we 
need to calculate this term if the pressure coefficient is required on each of the upper 
and lower surfaces. To calculate the local rate of change of the disturbance velocity 
potential, we need to know the velocity potential at this location at two successive time 
steps by integrating the velocity from an undisturbed position. Because of the long 
computational time involved, we only calculate the net pressure force on the wing in the 
unsteady-flow problem. However, calculations of the pressure distribution on the upper 
and lower surfaces are now under consideration because it is necessary for any boundary- 
layer calculations. 


The net pressure coefficient is given by 

• &t^.) X f. - 




(19) 


The total-load coefficients are obtained by integrating the net pressure coefficient on 
the wing. 
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V. NUMERICAL RESULTS 


A computer code which accounts for the g6’er5»l fcrnrilafion presented here is deve- 
loped. The code can be used to solve steady, unsteady, symmetric or asymmetric flow 
problems. Typical results obtained with this code are presented in this section. The 
required computation time and convergence studies are discussed in references 7-11 and 
13-16. 

1 . Steadily Rolling Wing at Zero Angle of Attac k 

In Figures 3-5 we show typical solutions of the free-vortex lines of a delta wing 
(AR = 0.7) at different rolling velocities 0.2, 0.4 and 0.6. In each figure, a three- 
dimensional view and a plan view are given. The traces of the helical vortex cones are 
indicated on the upper view which show that the cones are antisymmetric. The size of the 
cone increases and it rolls-up more rapidly with increasing the rolling velocity. This 
is in agreement with the experimental results^ 

Figure T shows the spanwise pressure distribution on the upper and lower surfaces at 
the chordv'isu s’:-Llon x/c^ = 0.778, the present results are compared with the experi- 
mental resui! r. Harvey®^ and the theoretical models of references IS and 20. A re- 
markable agreement can be seen between the present results and the experimental ones. 

The methods of references 18 and 20 obviously overestimate the suction pressure peak and 
this is attributed to the simplified representation of the separated flow. 

In Figures 7-12 v/e chow the effect of the yaw angle on the free-vortex lines and 
pressure distributions of a delta wing (AR = 0.7, a = 15°). 

The plan views in Figures 7, 9 and 11 show that the size of vortex cone emanating 
from the windward side of the wing increases while the size of the cone from the leeward 
side decreases. The former cone moves inboard while the latter cone moves outboard. 

The corresponding effect on the spanwise pressure distribution at the chordwise 
station x/c^ = 0.395 and comparisons with the results of Pull in*’ and the experimental 

results of Harvey” are given in Figures 8, 10 and 12. When the angle of yaw increases, 
the suction pressure peak on the windward side increases while that on the leeward side 
decreases. The present results are in a good agreement with the experimental ones. 

Pull in's results overestimates the suction peak on the windward side of the wing and it 
shows that the error increases with increasing the angle of yaw. 

3 . Unsteadily _Rol J ing W i ng at Zero Angle of Attack 

In Figure 13 we show the free-vortex lines at three successive time steps for a 

delta wing undergoing an unsteady rolling motion given by to = 0.6 + 0.1 sin(iTt/6). The 

X 

starting vortices shed from the edges of separation and later convected downstream can be 
seen in each view. 

Figure 14 shows the spanwise pressure distribution at the chordwise state x/c^ = 

0.65 at different time steps. The variation of the rolling-moment coefficient C with 

A 

the frequency of the sinusoidal variation of the rolling velocity Q is given in Figure 
15. In the initial time steps, we notice that the rolling-moment coefficient decreases 
as the frequency increases. Finally, the variation of the rolling-moment coefficient 
with the rate of roll is shown in Figure 16. 


13 



V I . CO NCLUDING REMARKS 

We have presented a general, nonlinear, discrete-vortex method for the lifting 
surfaces. The method is applied to delta wings undergoing different motions. The 
results are compared with several results of other investigators and with the experi- 
mental results of Harvey. They show that the present method provides accurate results 
which could not be obtained by the existing approximate methods. Moreover, the present 
method is not restricted by any geometrical parameter of the lifting surface. Moreover, 
it is exact because it is free from any small-disturbance assumption or any corresponding 
1 inearization. 
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b. Discrete-vortex system 
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Fiq. 2. Construction of the discrete-vortex system. 
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Fig. 11. Wake shape of a yawed delta wing in a steady flow, 
a = 15“, B = 10°, 8x8 bound lattice, AR = 0.7. 
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Fig. 14. Variation of the net spanwise pressure distribution 
with time on an unsteadily rolling delta wing, 

01 = - 0.4 - 0,1 sin(’irt/4), a = 0°, 6x6 bound 

A 

lattice, AR = 0.7, x/c„ = 0.65. 
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Fiq. 15. Variation of the rolling-moment coefficient with the 

frequency of rolling velocity of an unsteadily rolling 
delta wing, w = - 0.4 - 0.1 sinQt, a ~ 0°, 6x6 

A 

bound lattice, AR - 0.7. 



/[l/2yO U^S{b/2)] 









1. Report No. 2. Oovernment Accession No. 

NASA CR-1 45235 

3. Recipient's Catalog No. 

n. Title and Subtitle 

THREE DIMENSIONAL STEADY AND UNSTEADY ASYMMETRIC FLOW 
PAST WINGS OF ARBITRARY PLANFORMS 

5. Report Date 

September 1977 

6. Performing Organliation Code 

7 Authorfsl 

Kandil, 0. A.; Atta, E. H.; and Nayfeh, A. H- 

a. Performing Organiiation Report No 

10 Work Unit No 

505-15-22-01 

9 Purfurming O^Jan^^atlon Name and Address 

Virginia Polytechnic Institute and State University 
Blacksburg, Virginia 24061 

11 Contract or Grant No, 

NGR 47-004-090 

13 Ty[;e of Report and Period Covered 

Contractor Report 

12- Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

14 Sponsoring Agency Code 

1 

1 15. Supptememary Notes 

Grant Technical Monitor, Dr. E. Carson Yates, Jr. 

1 

16. Abstract j 


The nonlinear discrete-vortex method has been extended to treat the problem of 
asymmetric flows past a wing with leading-edge separation, Including steady and 
unsteady flows. The problem is formulated in terms of a body-fixed frame of reference, 
and the nonlinear discrete-vortex method is modified accordingly. Although the method 
is general, only examples of flows past delta wings are presented due to the avail- 
ability of experimental data as well as approximate theories. Comparison of our 
results with the experimental results of Harvey for a delta wing undergoing a steady 
rolling motion at zero angle of attack demonstrates the superiority of the present 
method over existing approximate theories in obtaining highly accurate loads. 

Numerical results for yawed wings at large angles of attack are also presented, jn 
all cases, total-load coefficients, pressure distributions and shapes of the free- 
vortex sheets are shown. 


17 Key Words ISuqgesied by Auihorlsl' 

Unsteady Aerodynamics 
Nonlinear Aerodynamics 
Vortex Flow Separation 
Leading-Edge Flow Separation 

Ifl Oisrrib'jtion 

Unclassified - Unlimited 

19 SecuYfty Classif fof t^irs rcportl ] 20 Security CMssif lof p.iqel 

Unclassified i Unclassified 

^ 

21 No. of Pages 

32 

22. Price’ 

$4.00 


* For sale by the National Technical information Service. Springfield Vngir-ia 22)6! 




















